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Abstract 

By transforming an exact stability condition, a new Nyquist-like plot is proposed to predict occur- 
rences of three typical instabilities in DC-DC converters. The three instabilities are saddle-node bifurcation 
(coexistence of multiple solutions), period-doubling bifurcation (subharmonic oscillation), and Neimark 
bifurcation (quasi-periodic oscillation). In a single plot, it accurately predicts whether an instability occurs 
and what type the instability is. The plot is equivalent to the Nyquist plot, and it is a useful design tool 
to avoid these instabilities. Nine examples are used to illustrate the accuracy of this new plot to predict 
instabilities in the buck or boost converter with fixed or variable switching frequency. 
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I. Introduction 

Operations of nonlinear DC-DC converters in continuous conduction mode (CCM) can be described 
exactly by the switching system shown in Fig. [T] (notations explained later) |[l], which leads to a sampled- 
data dynamics O]. Instability occurs when there exists a sampled-data pole outside the unit circle in 
the complex plane. There are three ways that the sampled-data pole leaves the unit circle, thus causing 
three typical instabilities in DC-DC converters |[T], ||2], iQ. When the sampled-data pole leaves the unit 
circle through 1 in the complex plane, the instability is generally a saddle-node bifurcation (SNB) ||2l 
(or pitchfork and transcritical bifurcations |@] which are less seen in DC-DC converters). The SNB 
generally involves coexistence of multiple solutions ||5], or sudden disappearances or jumps of steady- 
state solutions When the pole leaves the unit circle through -1, the instability is a period-doubling 
bifurcation (PDB) 0, which generally imolwes fast-scale subharmonic oscillation. When the pole leaves 
through a point other than 1 or -1 on the unit circle, the instability is a Neimark-Sacker bifurcation (NSB), 
which generally involves a slow-scale quasi-periodic oscillation lU, 13, ||3l. Many instability examples 
of DC-DC converters are shown in lITl. Il8l. Il9l. lITOl. ifTTI. 

Averaged models are traditionally applied to analyze DC-DC converters. It has been known that most 
averaged models cannot accurately predict the occurrence of PDB (subharmonic oscillation) By 
considering the sampling effects and increasing the system dimension, improved averaged models can 
predict the occurrence of PDB |[T2| . |[T3l . However, there is no single design-oriented tool like the Nyquist 
plot to predict the occurrences of all three instabilities. In a recent paper [|14i , two different approaches 
for slow-scale instability (NSB) and fast-scale instability (PDB) respectively are combined to predict 
instabilities in DC-DC converters. 

The switching system of Fig.[T]is a well-accepted model to describe the operation of DC-DC converters. 
The stability of the switching system of Fig. [T] is a pure mathematical problem and has a closed-form 
solution. It is well known |[T5l that the condition of pole locations can be easily transformed to a 
characteristic equation, and the Nyquist plot determines the stability. Then, one has the next question: 
how to derive an exact dynamics to be used for the Nyquist plot? 

In this paper, a new Nyquist-like plot, designated as "F-plot", based on sampled-data dynamics is 
proposed to predict the occurrences of all three instabilities in a single plot. The F-plot is based on the 
exact switching system. Therefore, the F-plot is also exact and it predicts the exact instability boundaries 
of these instabilities. It is applicable to converters under either current mode control (CMC) or voltage 
mode control (VMC). 

The remainder of the paper is organized as follows. In Sections |ll] and |llll the exact switching system 
and the sampled-data analysis are presented. In Sections JV] and |Vl a new F-plot is proposed and it 
is equivalent to the Nyquist plot. In Sections rvHlVIII[ the F-plot or Nyquist plot is applied to various 
examples to predict PDB, SNB, and NSB, respectively. In Section JXl similar methodology is applied to 
constant on-time control (COTC). Conclusions are collected in Section |Xl 

II. Brief Review of Exact Switching System 

A unique aspect about the switching system of Fig. [T] is that the ramp in VMC or CMC is represented 
by the same signal h(t) HI, HH, ^j}, GH. The operation of VMC and CMC is briefly reviewed here to 
make this paper self-contained. In VMC, the reference signal Vr controls the output voltage Vq- In CMC, Vr 
is generally replaced by v to denote a control signal to control the (peak) inductor current i^. Denote the 
source voltage as Vs. In the model, ^i, A2 G R^^^, Bi, B2 G R^''^ C, £^1, -E2 G R^^^, and D G R^^^^ 
are constant matrices, where N is the system dimension. Within a clock period T, the dynamics is switched 
between two stages, and 82- Switching occurs when the ramp signal h{t) := Vi + {Vh — mod T) 
intersects with the compensator output y := Cx + Du G R, where h{t) varies from a low value to 
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Figure 1. Exact switching system for a DC-DC converter in CCM. 



a high value Vh. Denote the ramp amplitude as Vm = Vh — Vi, and let the ramp slope be h{d) = ma 
Denote the switching frequency as fs := 1/T and let uJs '■= 27r/s. 



III. Small-Signal Analysis 

Let the steady-state duty cycle be D and let d := DT. Let the T-periodic orbit be x^{t). Let x^{d~) = 
Aix^{d) + Biu and x^{d'^) = A2X^{d) + B2U denote the time derivative of x^{t) at t = d~ and d'^, 
respectively. Let y^{t) = Cx^{t) + Du. One has y^{t) = Cx^{t). Confusion of notations for capacitance 
C and duty cycle D with the matrices C and D can be avoided from the context. 

In steady state, 



lU (1) 

T-d 



Jo 

fT-d 

xO(0) = e^^(^-^)x°((i) + / e^"'daB2U (2) 
y^{d) = Cx°{d)+Du = h{d) (3) 

Solving ([T])-© gives D and 

Generally the controller includes an integrator (with a pole at zero), making Ai (or A2) and I — e^^'^ 
non-invertible. In that case, the pole at zero can be replaced by a very small number, then Ai and /— e^^^ 
are invertible. Therefore, the invertibility of ^1 or / — e^^^ is not critical and can be resolved. This 
statement about invertibility of a matrix is not repeated later. 

Using a hat " to denote small perturbations (e.g., Xn = Xn — x'^{0)), where Xn is the sampled state at 
t = nT. From [H, |[T6l . |[T9l . ||20l . the linearized closed-loop dynamics is 

^ ^ y^{d-)-h{d) ' 

Instability occurs when there exists an eigenvalue z (also the sampled-data pole) of $ outside the unit 
circle. 
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Figure 2. A plant with unity negative feedback. 



IV. Discrete-Time Nyquist Plot 

Let $ = $0 - r^, where 

r = e^'^^-'^\x°{d-)-x°{d+)) (6) 
^ = (7) 

The dynamics (H)) can be transfomied into a plant (having an input dn and an output Wn, shown in Fig. |2]) 
with a unity negative feedback (dn = —Wn) 

Xn+l = ^QXn + Tdn 

(") 

dn = 

Although it looks like that ^ is arbitrarily transformed from dljl, the dynamics dU) does have a meaning. 
It is exactly the power stage (but including the compensator) dynamics |[20l with a PWM modulator gain 
l/{Cx^{d~) — ma), or l/{Cx^{d~) — ma)T if the duty cycle is used as the input. The modulator gain 
agrees with |[T2| which can predict PDB. It differs from l/niaT based on traditional averaging ll2n . ll22l 
which cannot predict PDB. 

The past research generally focuses on averaged dynamics. The Nyquist plot of the averaged dynamics 
without considering the sampling effect cannot predict PDB. Also, the past research focuses on dU, 
without transforming it into dH) which could lead to the Nyquist plot. 

Let the loop gain transfer function of dH) be J\f{z). The characteristic equation is 

M{z) = ^{zl -^qY^T = -I (9) 

The (discrete-time) Nyquist plot is J\f{e^^'^) as a function of u swept from to oJs/'^ for a half plot (or 
from —ijJs/2 to Ug/I for a full plot). 



V. Exact Stability Condition and the "F-Plot" 

Although using the Nyquist plot is enough to predict instability, a new plot is proposed to give a 
different perspective. Based on [[T] p. 46] and ll23l . lITSl . suppose z is not an eigenvalue of $o, then z is 
an eigenvalue of $ if and only if the following critical condition holds: 

M{z) := Cx^{d~) + Ce^i'^(z/ - $o)"^r = (10) 
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The proof is as follows. Suppose z is not an eigenvalue of $o- Then 

det[z/-$] = det[z/-$o]det[/+ (z/- «>o)~^r^'] 
= det[z/- $o](l + ^(^/- «'o)~^r) 

Then det[z/ - = requires that N(z) = ^(z/ - ^>o)~^r = -1, which leads to ^ by using ©-dTJ. 
Note that (ITOl ) is equivalent to the characteristic equation (|9]l through a scaled translation: 

M{z) = (Ci^id^) -ruaWiz) + Cx°{d^) (11) 

Although the stability condition (fTOl) looks unfamiliar, a special case of (fTOl) for PDB in CMC would be 
familiar. When the switching frequency is high, e"^^^ w gA2(T-d) ^ j -pjigj^ (-j^g critical condition ([TOl l 
with z = — 1 (for PDB) becomes 

^(CxO(d-) + Cx°(d+)) = + y°(d+)) = (12) 

In CMC, the (inductor current) slopes y^{d~) and y^{d~^) are generally expressed in most textbooks 
|[2T1 p. 448] as —mi and m2, respectively. Then, (fT2]) corresponds exactly to the well known minimum 
ramp slope ma = (m2 — mi)/2, or equivalently (ma — 1112) /{nia + mi) > — 1, required to stabilize the 
converter ll2n . ll24l . From (@), ([TOl i. or (fTZi . one sees that the stability is closely related to the signal 
slopes im, II20I. 

Instability occurs when the pole z leaves the unit circle, or equivalently, z satisfies the critical condition 
( fTOl ). One can define a plot based on the critical condition dTOl ) with z being swept around the unit circle. 
Let z = e*^ (on the unit circle), and define an "F-plot" in the complex plane as 

F{e) = Mie''^) = {CiP{d-) - niaWie'^) + C±^{d~) (13) 

which is a scaled translation of Nyquist plot. From ([T3T l. if A/'(e*^) = —1, one has F{6) = ma- The "-1 
point" in the Nyquist plot is translated to the "m^ point," (ma,0), in the F-plot. 

The idea behind the F-plot is exactly the same as Nyquist plot. As z is swept along the unit circle (and 
6 := ujT is swept from — vr to vr), one can determine the number of poles outside the unit circle (and 
hence the stability) based on the F-plot. Like the Nyquist plot, Im[F(0)] = —Im[F{—9)] and the F-plot 
is symmetric with respect to the real axis. One can make the F-plot for 9 G [0,7r] instead of [— 7r,7r]. 

Let P be the number of open-loop poles of the plant in Fig. |2] outside the unit circle. Generally, all 
eigenvalues of $0 are inside (or on) the unit circle 11201 . then P = 0. Let Z be the number of roots 
of the characteristic equation M{z). Let N be the number of (clockwise) encirclements of Nyquist plot 
A/'(e-''^^) around the -1 point. From |[T5l . N = Z — P = Z, and the following numbers are also equal to 
N: 

1) The number of poles for the dynamics dU outside the unit circle. 

2) The number of poles for the dynamics (|8]l outside the unit circle. 

3) The number of eigenvalues of <^ = <I>o — T^* outside the unit circle. 

4) The number of encirclements of the F-plot around the nia point. 

In terms of the pole z for @ or dSJ, the pole z satisfies the following equivalent equations 

1) det[zl - = det[z/ - $0 + T^*] = 0, 

2) the characteristic equation, M{z) = '^{zl — <l>o)~^r = —1, and 

3) the critical condition, M[z) = Cx^{d~) + Ce^^'^{zl - $o)~^r = m^. 

The critical conditions of PDB, SNB, and NSB are summarized in Table IH For example, PDB occurs 
when A = — 1, F{'k) = nia, or A/'(— 1) = — 1. The procedure of stability analysis is as follows. First, solve 
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Table I 

Critical conditions of PDB, SNB, and NSB. 
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Figure 3. Buck converter under PVMC. 



([T])-© and obtain D and x^{d), which can be also obtained approximately from the average analysis. 
Then, make the F-plot based on (fTOb and ([T3] ). Using ([S]!-©, one can also make Nyquist or Bode plot. 

In the following nine examples, various control schemes (such as average current mode control and 
type-III compensator) for buck or boost converters are used. The F-plot is compared with Nyquist and 
Bode plots. All of these plots accurately predict the three typical instabilities. 

The F-plot starts from F{0) (marked as *) and ends at F{tt) (marked as o). The other half of symmetric 
part for 6 G [— vr, 0] is not plotted for briefness. The nia point, located at (m^, 0) in the complex plane, 
is marked as +. In each example, two F-plots (as solid and dotted lines, respectively) are made before 
and after the bifurcation to show the onset of different instabilities predicted by the F-plot. 

VI. Prediction of Period-Doubling Bifurcation (PDB) 

Example 1. (The F-plot accurately predicts PDB and shows the required ramp slope to stabilize the 
converter.) Consider a widely studied buck converter with proportional voltage mode control (PVMC) 
from |I6] shown in Fig. |3] The proportional feedback gain is kp = 8.4. The converter parameters are 
T = 400 /xs, L = 20 mH, C = 47 pF, R = 22 n, Vi = 3.8, 14 = 8.2, and Vr = 11.3 V. The PDB is 
known to occur when Vg = 24.5. 

Time-domain simulation. First, let Vg = 24. The converter is stable with a T-periodic orbit shown 
in Fig. m^a). Next, let Vs = 25. The T-periodic orbit is unstable. Subharmonic oscillation occurs with a 
2r-periodic orbit shown in Fig. IHb). 

F-pIot. Here, the ramp slope is rua = Vm/T = 11000. First, let Vg = 24, the converter is stable. 
Solving ([Hi-® gives D = 0.5. The F-plot (Fig. [5]) does not encircle the rria point (if it is plotted for 
9 € [— vr, vr] which includes the other half of symmetric part of the F-plot). At Vg = 24.5, F{7r) touches 
the nia point and PDB occurs. 
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Next, let Vs = 25, the converter is unstable. Solving (H}-© gives D = 0.48. The F-plot (Fig. [5]) 
encircles the nia point (if it is plotted for 9 G [— 7r,7r]). The F-plot accurately predicts the occurrence 
of PDB. The F-plot also shows how to adjust the ramp slope to stabilize the converter with Vg = 25. 
The F-plot in Fig. [5] shows that if nia > F{'k) = 11276 the converter is stabilized. To keep the original 
operating condition at Vg = 25 and D = 0.48, let Vi = 3.6856 and Vh = 8.3056. Now, nia = 11550. 
The sampled-data poles are —0.8202 it 0.0803j (inside the unit circle), confirming the stability. 

Nyquist/Bode plot. First, let Vs = 24. The Nyquist plot (Fig. [6l) does not encircle the -1 point. The 
Bode plot (Fig. |7]) shows a very small gain margin of 0.0933 dB. 

Next, let Vs = 25. The Nyquist plot (Fig. (6]) encircles the -1 point. The Nyquist plot starting from 
AA(1) (for a; = 0) to ^f{-l) (for uj = uJs/2) is the lower half in Fig. |6l Here, 7V(-1) is closer to the 
-1 point than J\f{l), indicating that it is a PDB instead of an SNB. The Nyquist plot accurately predicts 
the occurrence of PDB. The Bode plot is similar to Fig. |7]but with a gain margin of -0.108 dB at the 
subharmonic frequency ujs/2 = 7854 rad/s. □ 



15 20 25 



(a) (b) 
Figure 6. (a) Nyquist plot, (b) blown-up view of (a); solid line for 113 — 24 (stable) and dotted line for Vs ~ 25 (unstable). 




Figure 7. Bode plot for Vs — 24. 



The next example shows that PDB occurs not only in peak current mode control but also in average 
current mode control (ACMC). The PDB instability in ACMC was reported ten years ago ||25l but it 
was unnoticed. The instability phenomenon is very interesting because the PDB instability occurs only 
in a window value of a pole. Since the value of the pole does not affect the duty cycle but affects the 
ripple size, this example also shows that the PDB instability is unrelated to the ripple size. For details 
see Q, m, ifTTl . The instability is verified by three independent approaches: time-domain simulation, 
the sampled-data pole trajectories, and F/Nyquist/Bode plots. 

Example 2. {The F-plot accurately predicts the subharmonic oscillation window in a buck converter 
under ACMC.) Consider an ACMC buck converter ||26] p. 114] shown in Fig. |8] The ACMC type-II 
compensator is 

KAl + — ) 



The converter parameters are Vs = 14 V, = 5 V, Vr = 0.5, fs = 50 kHz, L = 46.1 /xH, C = 380 
/uF with ESR Rc = 0.02 VL, R = I ^, = 0.1 = 0, V^=l, = 50000, = 75506, and 
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ujz = 5652.9 rad/s. Let (5 = 1 be the integrator pole which is small and it does not affect much the 
dynamics. 

Time-domain simulation. The compensator pole Up is varied from 0.14tJs to O.Slwg. An unstable 
window of ojp between O.lSa;^ and 0.49ti;s was found and reported in ||25| . When ujp is inside the window, 
the subharmonic oscillation occurs. 

For example, let ujp = 0.15a;s. The converter is stable (Fig. |9]l. Let ojp = 0.49a;s- The converter is 
unstable with subharmonic oscillation (Fig. [TOll. Next, let ojp = O.SlWs. The converter is stable again 

(Fig. m- 

Independent sampled-data analysis. The sampled-data pole trajectories for ojp/ujs £ (0.1,0.8) are 
shown in Fig. [121 There are four poles. Two poles are fixed around 0.88, and 0.95 (« e~Rc). A pole 
leaves the unit circle through -1 when ujp = O.lSa;^, and enters the unit circle when ujp = 0.49a;s- This 
explains exactly the instability window of ujp. 

F-plot. Solving O-© gives D = 0.3571. First, let bjp = O.lSt^s, the converter is stable. The F- 
plot is shown in Fig. [131 The non-smoothness of the curve is due to the small integrator pole 6. For 
6 € (0, vr], the curve is actually smooth. When 9 is close to 0, the F-plot trajectory suddenly jumps from 
F(0) = — 5 X 10^ to a large value and forms an abrupt bent. The blown-up views of the F-plot are shown 
in Fig. [I4l The F-plot does not encircle the nia point. The sampled-data poles are 0.8868, 0.9521, and 
—0.6406 lb 0.1557j (inside the unit circle), also confirming the stability. 

Next, let ujp = 0.490;^, the converter is unstable. The F-plot is shown in Fig. [151 and it encircles the 
nia point. The sampled-data poles are -1.0080, -0.0511, 0.8826, 0.9533, also confirming the instability. 

To show the instability window, let ojp = O.Sa;^, the converter is stable again. The F-plot is shown 
in Fig. [161 and it does not encircle the nia point. The sampled-data poles are 0.0487, 0.9936, 0.8825, 
0.9533, also confirming the stability. The F-plot accurately predicts the instability window. 

Nyquist/Bode plot. First, let Up = Q.lbuJs- The Nyquist plot (Fig. Wl\ does not encircle the -1 point. 
The Bode plot (Fig.[T8l) shows a positive gain margin. Second, let ojp = 0.49^1;^. The Nyquist plot (Fig. [T9l ) 
encircles the -1 point. The Bode plot (Fig. l20l ) shows a negative gain margin. Third, let = 0.5tJs- 
The Nyquist plot (Fig. |2TI ) does not encircle the -1 point. The Bode plot (Fig. l22l) shows a positive gain 
margin. □ 

Example 3. (With phase margin of 38.9° based on the averaged model, subharmonic oscillation still 
occurs.) A converter with a type-Ill compensator is shown in Fig. |23l A type-III compensator ll27l p. 




Figure 8. System diagram of an ACMC DC-DC converter. 
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Figure 9. Stable T-periodic solution, ojp = O.lSuja. 




Figure 10. Subharmonic oscillation, Up = 0.49a;s. 




Figure 11. Stable T-periodic solution, LOp = O.Slojs. 



261] has three poles, two zeros, and a gain Kc, with a transfer function 
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A typical guideline |[27l p. 412] popular in industry to set the parameters of the compensator is as 
follows. Set one pole at 6 ^ (as an integrator), and set p2 = 1/RcC. Set the gain Kc to adjust the 
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Figure 15. The blown-up view of F-plot sliows that the rUa point is encircled by the F-plot, cjp = 0.49tJs 
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Figure 16. The blown-up view of F-plot shows that the nia point is not encircled by the F-plot, ujp = 0.5a;s 




Figure 17. 
-1 point. 



(a) The Nyquist plot for ujp 



— 0.15u}s, (b) blown-up view of (a) shows that the Nyquist plot does not encircle the 
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Gm = 0.47 dB {at 1 .57B+005 rad/sec) , Pm = 5.48 deg {at 1 .33B+005 rad/sac| 

60 I ■ ■ ■ ■ — 1— ■ — — 




10^ 10'' lo"" 10^ 

Frequency (rad/sec) 



Figure 18. The Bode plot shows a positive gain phase margin, Up = O.lScJs. 




(a) (b) 
Figure 19. (a) The Nyquist plot for ujp — O.lScJs, (b) blown-up view of (a) shows that the Nyquist plot encircles the -1 point. 




Figure 20. The Bode plot shows a negative gain phase margin, ujp = 0.49a;3. 



phase margin and the crossover frequency. Let zi = Kz/vLC and Z2 = 1/vLC, where is a zero 
scale factor to have additional flexibility to adjust the phase margin and the crossover frequency. The 
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(a) (b) 

Figure 21. (a) The Nyquist plot for oip — O.lSa;^, (b) blown-up view of (a) shows that the Nyquist plot does not encircle the 
-1 point. 




Figure 22. The Bode plot shows a positive gain phase margin, Up = O.So^s. 



zero scale factor used in industry typically varies between 0.1 and 1.2. As will be shown later, a 
smaller value of may lead to subharmonic oscillation. Taking into account the above guidelines, the 
compensator has a transfer function 



Consider a buck converter with the type-Ill compensator (1161) . Exactly the same parameters as in a 
technical note EH are used: = 1/T = 300 kHz, L = 900 nH, C = 990 ^F, R = 0.4 Q, = 5 mQ, 
Vr = 3.3 V, Vm = 1.5 V, Kc = 7.78 x 10^, zi = 1/2/LC = 1.675 x 10^ Z2 = l/y/IC = 3.35 x 10^, 
Up = Us/2 = 9.425 X 10^, and p2 = 1/RcC = 2.02 x 10^. Let 5 = 1 be the integrator pole which is 
small and it does not affect much the dynamics. 

Simulation (Fig. l24b shows that subharmonic oscillation occurs when = 16 V (Z) « 0.206). This is 
also confirmed by the exact sampled-data analysis with a sampled-data pole at - 1 when the subharmonic 
oscillation occurs. Based on the averaged model for Vs = 16, the loop gain frequency response (Fig. 



15 




Figure 23. A DC-DC converter with a type-III compensator Gc(s) = Hn{s) / Hd{s). 




Figure 24. Signal waveforms showing the subharmonic oscillation. 



shows a phase margin of 38.9°. The frequency response also shows an infinite gain margin because the 
phase never reaches -180°, which means no matter how much Vg increases (to increase the loop gain), 
the converter is expected to be stable based on the averaged model. However, subharmonic oscillation 
still occurs when = 16. □ 

Example 4. {Unstable window of pi, unrelated to the ripple size of y^{t).) Consider again Example |3] 
where pi = ixj^ll- Now, vary p\ from ^Xug to CGw^. 

Unstable window. An unstable window of p\ G (0.23, 0.5)^;^ is found. The value of p\ adjusts the 
ripple size of ]p{t\ A larger p\ leads to a larger ripple of ipii). In |[T4l . it is hypothesized that the ripple 
size of ipii) is related to subharmonic oscillation. The following simulation shows that the ripple size 
of (t) is unrelated to subharmonic oscillation. 

Time-domain simulation. The unstable window of pi is confirmed by time-domain simulation. For 
pi = 0.2a;s, the ripple size of y is small, and the converter is stable (Fig. |26]l. For pi = Q.2Auis, the 
ripple size of y is larger, and the converter is unstable (Fig. ^7). For pi = Q.Qujs, the ripple size of y 
is even larger, but the converter is stable (Fig. |28}. Comparing Figs. |26ll28l the ripple size of y^{t) is 
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Figure 25. Loop gain frequency response shows a phase margin of 38.9° and an infinite gain margin based on the averaged 
model, but the subharmonic oscillation still occurs. 




Figure 26. The converter is stable, pi = Q.2us. 



unrelated to subharmonic oscillation. This shows a counter-example for the hypothesis proposed in ||T41 
that the ripple size of y'^(t) is related to the subharmonic oscillation. 

Confirmation by sampled-data pole trajectories. The unstable window of pi is also confirmed by 
the sampled-data pole trajectories. The sampled-data pole trajectories for O.loJs < pi < CGw^ are shown 
in Fig. |29l Three poles are fixed around 0.9485, 0.8853, and 0.51. A pole leaves the unit circle through 
-1 when pi = 0.23a;s, and enters the unit circle when pi = Q.buJs- This explains exactly the unstable 
window of pi. 

Confirmation by F/Nyquist/Bode plots. Here, nia = Vm/T = 45000. First, let ojp = Q.2uja, and the 
converter is stable. The F-plot (Fig. |30] l does not encircle the nia point. The Nyquist plot (Fig. [3T]) does 
not encircle the —1 point. The Nyquist path goes to infinity and connects to form a loop (not shown). 
This statement is not repeated later when the Nyquist path goes to infinity. The Bode plot (Fig. [32] ) shows 
positive gain/phase margins. 

Second, let ojp = 0.5ws, and the converter is unstable due to PDB. The F-plot (Fig. [33]) encircles the 
nia point. The Nyquist plot (Fig. [34] ) encircles the —1 point. The Bode plot (Fig. [35] ) shows a negative 
gain margin. 
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Figure 27. Subharmonic oscillation, pi = 0.24aJs. 
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Figure 28. The converter is stable with a larger ripple of y^{t), pi = 0.6ujs. 
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Figure 29. Sampled-data pole trajectories for O.Iljs < pi < 0.6a;s. 

Third, let Up = O.Gws, and the converter is stable. The F-plot (Fig. |36]l does not encircle the ma 
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(a) (b) 

Figure 30. (a) F-plot for 6 G [0, tt], uip = Q.2u)s\ (b) blown-up view of (a) shows that the F-plot does not encircle the m,a 
point. 





(a) (b) 
Figure 31. (a) Nyquist plot, u}p = 0.2tJa; (b) blown-up view of (a) shows that the Nyquist plot does not encircle the -1 point. 



point. The Nyquist plot (Fig. [37l l does not encircle the —1 point. The Bode plot (Fig. [38l) shows positive 
gain/phase margins. 

In this example, the unstable window of pi is verified by three different approaches: time-domain 
simulation, the sampled-data pole trajectories, and F/Nyquist/Bode plots. □ 



VII. Prediction of Saddle-Node Bifurcation (SNB) 



Example 5. {The F-plot accurately predicts SNB.) Consider a boost converter from |^9l . ll30l shown 
in Fig. |39] The converter parameters are T = 2 /iS, = 4 V, L = 5.24 ^uH, C = 0.2 /iF, 7? = 16 0, 
Vt = 0.48, = 0, and Vh = 1. Let the feedback gains for the current and voltage loops respectively 
be ki = —0.1 and = 0.01. The switch is on when h{t) < y := Vr — hii — k^Vo and it is off when 
otherwise. 
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Figure 32. The Bode plot shows positive gain/phase margins, ujp = 0.2ujs 
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Figure 33. (a) F-plot for 9 G [0, tt], cop = O.Scjs; (b) blown-up view of (a) shows that the F-plot encircles the rua point. 





(a) 



(b) 



Figure 34. (a) Nyquist plot, ujp = O.ScJj,; (b) blown-up view of (a) shows that the Nyquist plot encircles the -1 point. 



Bode Diagram 
Gm = -0.00053 dB (at 9.42e+005 rad/sec] , Pm ^ Inf 




Figure 37. Nyquist plot does not encircle the -1 point, ujp = 0.6cjs. 
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Figure 38. The Bode plot shows positive gain/phase margins, ujp = O.Gujs 
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Figure 39. Boost converter under multi-loop state feedback. 



Time-domain simulation. Three solutions coexist: a stable T-periodic orbit with D = 0.586 (solved 
from ©-©), an unstable T-periodic orbit with D = 0.71, and a stable DC solution with {iL,Vo) = 
{vs/tl,^), where is inductor resistance. The converter is thus bistable. The two T-periodic orbits are 
shown in Fig. |40l The final state is sensitive to the initial condition. First, let the initial condition be 
(iL(0), 1^0(0)) = (2.2, 16.2), close to the unstable T-periodic orbit shown in Fig. gO] In about SOT, the 
state trajectory moves toward the stable T-periodic orbit as shown in Fig. |4T] The stable orbit attracts the 
state trajectory. Next, let (i/,(0), i;o(0)) = (2.3, 16.2), still close to the unstable orbit. The state trajectory 
moves away from the unstable orbit as shown in Fig. Hi] agreed with BOL The unstable orbit expels the 
state trajectory toward the stable DC solution. After t > 8T, the switch is always on (since y{t) > h{t)), 
and the inductor current increases toward Vs/ri. 

The bifurcation diagram using Vr as the bifurcation parameter is shown in Fig. |42]with the DC solution 
omitted. One sees that SNB occurs at Vr = 0.496 and D = 0.65 when the two T-periodic orbits coalesce. 
For Vr > 0.496, no T-periodic orbit exists. For Vr = 0.48, the bifurcation diagram shows exactly the two 
T-periodic orbits with D = 0.586 and D = 0.71. 

F-plot. Here, the ramp slope is rua = Vm/T = 500000. First, let D = 0.586. The F-plot (Fig. |43l) 
does not encircle the ttIq point. The sampled-data poles are 0.8045 it 0.4510j (inside the unit circle), 
also confirming the stability. 
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Figure 40. Two coexisting T-periodic orbits in state space. 
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(a) (b) 

Figure 41. Stability is sensitive to the initial condition, (a) The state trajectory starting from {iL{0),Vo{0)) = (2.2,16.2) 
moves toward the stable orbit; (b) The state trajectory starting from (11,(0), Wo(0)) — (2.3, 16.2) moves away from the unstable 
orbit. 



Note that the F-plot (such as Fig. [43] ) associated with SNB is a little different from the F-plot (such 
as Fig. Is) associated with PDB. In Fig. [43l F{0), which is associated with SNB, is closer to the ma 
point than F{it), and the F-plot for 9 G [0, vr] is almost below the real axis. In Fig. [5j F(7r), which is 
associated with PDB, is closer to the lUa point than F(0), and the F-plot for 9 G [0,7r] is above the real 
axis. However, in the both plots, they are plotted clockwise as 9 increases. 

Next, let D = 0.71. The F-plot (Fig. [43l) encircles the nia point. The F-plot accurately predicts the 
occurrence of SNB. The sampled-data poles are 1.5891 and 0.6501, also confirming the instability. 

Nyquist/Bode plot. First, let D = 0.586. The Nyquist plot (Fig. |44(a)| ) does not encircle the -1 point. 
The Bode plot (Fig. |44(b)[ ) starts from -6 dB, and it increases to 2 dB with the phase decreased from 
180° (away from the -1 point, thus implying stability). 

Next, let D = 0.71. The Nyquist plot (Fig. |44(a)| ) encircles the -1 point. Different from Example 1, the 
Nyquist plot starting from M{1) (for lo = 0) to M{—1) (for io = u]s/2) is the upper half in Fig. |44(a) 



Here, M{1) is closer to the -1 point than J\f{—1), indicating that it is an SNB instead of a PDB. The 
Nyquist plot accurately predicts the occurrence of SNB. The Bode plot ((Fig. |44(b)| )) starts from 5 dB 
with the phase decreased from 180° (encircling the -1 point, thus implying instability). □ 
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(a) 

Figure 44. (a) Nyquist plot, (b) Bode plot, solid line for D 



(b) 

= 0.586 (stable) and dotted line for D = 0.71 (unstable). 
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Figure 45. A buck converter under current-mode control with CCL 



The next example from lITSl shows that using a constant current load (CCL, such as an LED) in 
a CMC buck converter may cause an SNB instability. The instability is verified by four independent 
approaches: bifurcation diagram, coexistence of multiple solutions, the sampled-data pole trajectories, 
and F/Nyquist/Bode plots. 

Example 6. (Current mode control (CMC) buck converter with CCL.) Consider a CMC buck converter 
(Fig. IS) with the following parameters IH]: = 1 A, = 10 V, = 1 MHz, Rc = 50 mO, L = 10 
fiR, and C = 20 fjF. No ramp is used here and nia = 0. 

Bifurcation diagram. Let ic (peak inductor current control signal) be the bifurcation parameter. The 
bifurcation diagram is shown in Fig. SNB occurs when v = 1.125, D = 0.4998, and Vq = 5. Note 
that the focus of this paper is on the T-periodic solution. Those non-T-periodic attractors are not shown 
in the bifurcation diagram to prevent detraction of the focus on SNB. The PDB occurs when D = 0.5006. 
As reported in fTSl, both PDB and SNB occur simultaneously around D = 0.5 if Rc = 0- 

Coexisting multiple solutions. Generally in SNB, there is a hysteretic loop as shown in Fig. |46] In 
the figure, the upper solid line is for the operation when the switch is always on (hence, D = 1), and 
the dashed line and the lower solid line are for unstable and stable T-periodic solutions respectively 
with duty cycle less than 1. For v < 1, the switch is always off (D = 0). For 1 < ic < 1.125, there 
are three solutions: one stable T-periodic solution (with D < 0.4998), one unstable T-periodic solution 
(with D > 0.4998), and the third (stable) DC solution being that the switch is always on (D = 1). Take 
ic = 1.12, for example. The two T-periodic solutions, with duty cycles 0.4 and 0.6, respectively, are 
shown in Fig. |47] The two solutions have the same peak inductor current 1.12, the same average inductor 
current and the same current ripple amplitude. 

Sampled-data pole trajectories. Each solution has two poles shown in Figs. |48] and |49l respectively. 
The stable solution (with D < 0.4998) has two sampled-data poles inside the unit circle. The unstable 
solution (with D > 0.4998) has two poles outside the unit circle for D > 0.5006 (where PDB occurs), 
and it has one pole outside the unit circle and one pole inside the unit circle for 0.4998 < D < 0.5006. 
When the converter operates with a periodic solution and v is increased a little above 1.125, the output 
voltage will jump up from 5 V to 10 V. Similarly, when the converter operates with D = 1 and ic is 
decreased a little below 1, the output voltage will jump down from 10 V to V. The jumping up and 
down forms a hysteretic loop shown in Fig. |46] 

If SNB is not noticed, one may believe the converter is stable for v < 1.125, which is wrong. For 
1 < ic < 1.125, the converter is actually bistable, which is a type of instability. There are two stable 
modes, one mode with D = 1 and the other with D < 0.5. The stable mode with D < 0.5 is T-periodic 
and is locally stable. One mode may jump to the other mode. This example illustrates the importance to 
uncover the existence of SNB. 
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Figure 46. Bifurcation diagram shows the occurrence of SNB at ic = 1.125, solid line for stable solution, dashed line for 
unstable solution 




Figure 47. Coexisting two steady-state inductor currents, solid line for stable solution, dashed line for unstable solution 
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Figure 48. The locus of the first pole as D varies 
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Figure 49. The locus of the second pole as D varies 
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(a) (b) 
Figure 50. (a) F-plot for £ [0,7r], D — 0.4; (b) blown-up view of (a) shows that the F-plot does not encircle (0,0). 



F-pIot. First, let D = 0.4, the converter is stable. Since no ramp is used, nia = 0. The F-plot (Fig. |50l ) 
does not encircle the rua point. 

Next, let D = 0.6, the converter is unstable. The F-plot (Fig. |5T) encircles the nia point twice (if a 
full plot for 6 G [— TT, tt] is made). The two encirclements are due to the fact that one pole is greater than 
1 (for SNB) and the other less than -1 (for PDB). Note that here F(0) ^ F{tt), and the double-loop 
F-plot encircles the point simultaneously, implying the simultaneous occurrence of SNB and PDB 
around D = 0.5. If F{0) and F{tt) are apart, then two encirclements of the F-plot around the nia point 
generally imply the occurrence of NSB with a pair of complex poles crossing the unit circle. 

Nyquist/Bode plot. First, let D = 0.4, the Nyquist plot (Fig. |52ll does not encircle the —1 point. The 
Bode plot (Fig. [53l ) shows positive gain/phase margins. Next, let D = 0.6, the Nyquist plot (Fig. |52j 
actually having two loops close to each other) encircles the —1 point twice. The Bode plot (Fig. [53l) 
shows a negative gain margin of -1.93 dB. □ 
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Figure 53. The Bode plot, solid line for D — 0.4 (stable) and dotted line for D — 0.6 (unstable). 
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ramp h{t) 




Figure 54. Buck converter with input filter. 




Figure 55. Stable T-periodic orbit, Rp = 34. 



VIII. Prediction of Neimark-Sacker Bifurcation (NSB) 

Example 7. {The F-plot accurately predicts NSB and the associated oscillation frequency.) Consider a 
buck converter with input filter from |l2] shown in Fig. |54]with Vs = 15.8 V. The power stage parameters 
are the same as in Example 1. The input filter parameters are Lf =2.5 mH and C/=160 /iF, and Rp is 
varied from 1 $7 to 100 Q. The NSB occurs when Rp = 38.85 O. 

Time-domain simulation. First, let Rp = 34. The converter is stable with a stable T-periodic orbit 
shown in Fig. [55] Next, let Rp = 39. The T-periodic is unstable. The state trajectory moves toward a 
quasi-periodic orbit shown in Fig. |56] The "period" is around lOT. The NSB may lead to a route to 
border-collision bifurcation or chaos IIBTI . The eigenvalue trajectory of <I> as Rp varies from to 100 is 
shown in 121 Fig. 5.43]. One pair of eigenvalues is almost fixed at — 0.5963±0.5301j, while the other pair 
moves as Rp varies. The NSB occurs around Rp = 38.85, when a pair of eigenvalues 0.8087 ± 0.5883j 
crosses the unit circle. 

F-plot. First, let Rp = 34. Solving O-© gives D = 0.79. The F-plot (Fig. [ST]) does not encircle the 
nia point (because if a horizontal line is drawn from the nia point to the right, it does not intersect with 
the F-plot). As Rp increases close to 38.85 when NSB occurs, F{9) = rua (and the F-plot touches the 
nia point), where 6 or it. 

Next, let Rp = 39. The F-plot (Fig. |57] ) encircles the nia point twice (because if a horizontal line 
is drawn from the nia point to the right, it intersects twice with the F-plot for 9 G [— 7r,7r]). The two 
encirclements agree with the fact that there are two unstable complex poles. The F-plot accurately predicts 
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the occurrence of NSB. 

The F-plot also predicts the oscillation frequency when the instability occurs. On the F-plot (Fig. [57] ). 
F(0.629) = 11010 — 464j is close to the nia point. The expected oscillation frequency, agreed with 
Fig. [56l is 0.629/r = 1573 ^ Ws/lO which is close to the resonance frequency of the input filter 
l/^/LjC]= 1581 rad/s. 

Nyquist/Bode plot. First, let Rp = 34. The Nyquist plot (Fig. |58]) does not encircle the - 1 point. The 
Bode plot (Fig. [59l ) shows positive gain and phase margins. Next, let Rp = 39. The Nyquist plot (Fig. [58l) 
encircles the -1 point twice. The Nyquist plot accurately predicts the occurrence of NSB. The Bode plot 
(Fig. [59l) shows a negative phase margin exactly at 1573 rad/s. The Bode plot intersects with the —180° 
line twice, also showing the two encirclements. □ 

Example 8. {Accurate prediction of NSB in a buck converter with a type-Ill compensator ) Consider 
again Example [3] where = 5 mU. Now, change it to R^, = 0.427 mf]. Based on the sampled-data 
analysis, the poles are -0.276 ± 0.9618j, 0.9477, 0.8884, and 0.0259. The complex poles outside the 
unit circle indicates the occurrence of NSB. The time-domain simulation (Fig. [60]) shows a quasi-periodic 
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(a) (b) 
Figure 58. (a) Nyquist plot, (b) blown-up view of (a); solid line for Rp = 34 (stable) and dotted line for Rp = 39 (unstable). 
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Figure 59. Bode plot (blown-up view), solid line for Rp — 34 (stable) and dotted line for Rp — 39 (unstable). 



solution with a "period" around 3.4T (corresponding to a frequency of 554400 rad/s), confirming the 
occurrence of NSB. The NSB may lead to a route to border-collision bifurcation or chaos lISTI . 

Confirmation by F/Nyquist/Bode plots. The F-plot (Fig. [6T]| encircles the rria point twice (if a full 
plot for 9 G [— vr,7r] is made). The Nyquist plot (Fig. [62] ) encircles the —1 point twice. The Bode plot 
(Fig. [63l) shows negative gain/phase margins. It also shows an oscillation frequency around 554400 rad/s, 
agreed with the time-domain simulation. The Bode plot intersects with the —180° line twice, also showing 
the two encirclements. □ 



IX. Similar Methodology Applied to Constant On-Time Control (COTC) 
For COTC, from |[TOl . the linearized sampled-data dynamics is 

Xn+l = ^Xn 



(17) 
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Figure 62. (a) Nyquist plot, Rc = 0.427 mQ; (b) blown-up view of (a) shows that the Nyquist plot encircles the -1 point. 
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Figure 63. The Bode plot shows negative gain/phase margins, Rc = 0.427 mQ. 



where 



dXn dTn dTn dXr. 



{x„,u„,Tr,)=ix°{0),u,T) 



Cx^{0~)-ma 



From ifTOl . one has the following result. Suppose that z is not an eigenvalue of e^^^'^ d)^Aid_ Then z 
is an eigenvalue of $ if and only if 

C{I - z-ie^^ V^(^-'^))-ii°(0-) = ma (19) 
The proof is as follows. Suppose z is not an eigenvalue of e^2(T-d)gyiid^ ^j^gj^ 
det[.2/-$] = det[z/-e^^(^-'^)e^^'^] • 

det[/ + {zl - eMT~d)^Mdyl, iy~)C yA.i^T-d)^Md 



(1 + Ce^-^-^^-'^^e^^\zl - ^^■'''-''^^'')-\c£^S^)^ 



Therefore, z is an eigenvalue of <^ if and only if 



1 + Ce^-^i'r-^)e^^\zl - ^A.i:r-d)^A.dyi> ■^"(O' J ^ ^ q 

Ci;'J(0 ) — ma 

which can be rearranged as 

Ci°(0-) + C(ze-^^(^-'^)e-^^'^ - /)-ix°(0-) = m,, (20) 

leading to (O based on the matrix equality / + {^ze~^'^^^"'^^ e~^^'^ - = ij - z~^e^i'^e^=(^~'^))^^ 
As discussed above, instability occurs when there exists an eigenvalue (also the sampled-data pole) z 
of ^ outside the unit circle. Let z = e*^ (on the unit circle), and denote the left side of ( fT9l ) as a function 
of 9, designated here as an "F-Plot" in the complex plane, 

F{e) = C{I - e-^^e^^ V^(^-'^))-ix°(0-) (21) 
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Figure 64. A voltage mode COTC buck converter. 



Like the Nyquist plot, Im[F(0)] = — Im[F(27r — 9)] and the F-plot is symmetric with respect to the 
real axis in the complex plane. One can make the F-plot for 9 G [0,7r] instead of [0,27r]. Also like the 
Nyquist plot, instability occurs when the F-plot encircles the point {nia, 0) {"rua point") in the complex 
plane. 

For stability analysis, the dynamics ([TTl i with ii„ = is equivalent to 

Un — 'if.Ln ■— Cx°{0-)-m^ 

or equivalent to an open-loop plant (with an input dn and an output "^Xn) as shown in (l22l ) with a unity 
negative feedback. Let the loop gain transfer function of (|22T i be J\f{z). One has M{z) = '^{zl — <^q)^^T . 
The dynamics (l22l) has a characteristic equation 1 + ^{z) = 0, which is equivalent to ( fT9l ). 

Example 9. {F/Nyquist/Bode plot accurately predicts PDB in COTC.) Consider a voltage mode COTC 
buck converter (Fig. |64]i from ll32l with the following parameters: t;,, = 5 V, T = 3 ^s, d = 1.2 ^s, 
R = 0.5 ^l, Rc = 20 mO, L = 2 fiU, and C = 20 ^F. Here, D = d/T = 0.4 and Vr = VsD = 2. Also 
no ramp is used here and nia = 0. With these parameters, instability occurs as shown in ||32] Fig. 9]. 

Sampled-data pole. Based on the exact sampled-data analysis by calculating the eigenvalues of 
one pole is and the second pole is -1.1, implying occurrence of PDB. 

F-plot, Nyquist plot and Bode plot. The F-plot is shown in Fig. |65] and it encircles the nia point (if 
it is plotted for 9 G [0, 2tt] to complete the other half of symmetric part). The F-plot accurately predicts 
the instability. The Nyquist plot is shown in Fig. [66] and it encircles the -1 point. The Bode plot is shown 
in Fig. |67] which shows a negative gain margin of -0.135 dB. □ 



X. Conclusion 

By transforming an exact stability condition, a new F-plot is proposed to predict occurrences of all 
three typical instabilities in DC-DC converters in a single plot. The F-plot predicts exactly the critical 
parameters (such as Vg = 24.5 in Example 1) when the instabilities occurs because it is based on the 
exact switching system of Fig. [T] The F-plot is equivalent to the Nyquist plot and therefore shares similar 
properties. It is symmetric with respect to the real axis in the complex plane and it shows instability if it 
encircles the nia point in the complex plane. The Nyquist plot or Bode plot of a particular discrete-time 
dynamics ([8]l, not reported before, also accurately predicts the instability. The F-plot offers a different 
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(a) (b) 

Figure 65. The F-plot starts from F(0) (marked as *) and ends at F{tt) (marked as o), (a) 9 £ [0, tt], (b) blown-up view of 
(a) shows that the nia point (marked as +) is encircled by the F-plot (if the lower symmetric half is plotted). 




(a) (b) 
Figure 66. (a) The Nyquist plot, (b) The blown-up view of Nyquist plot shows encirclement of the -1 point. 



perspective from the Nyquist plot or the Bode plot. Depending on the closeness of the nia point to either 
F(0), F{7r), or F{9) (where / or vr), one can predict respectively which instability (SNB, PDB, or 
NSB) is likely to occur. The F-plot also shows the required ramp slope to avoid these instabilities and 
predicts the oscillation frequency associated with NSB. Nine examples are used to show the accurate 
prediction of F/Nyquist/Bode plot to predict PDB/SNB/NSB in the buck or boost converter with fixed or 
variable switching frequency under various control schemes. 
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